Topological Force and Torque in Spin-Orbit Coupling System 
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The topological force and torque are investigated in the systems with spin-orbit coupling. It is 
demonstrated that the topological force and torque appears as a pure relativistic quantum effect 
in an electromagnetic field. The origin of both topological force and torque is the Zitterbewegung 
effect. Considering nonlinear behaviors of spin-orbit coupling, we address possible novel phenomena 
driven by the topological forces. 
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Recently, the spin-orbit coupling has become an in- 
teresting topic due to the spin Hall effect [1-5]. It pro- 
vides an efficient route to generate and control quantum 
spin state electrically. Generally speaking, the spin-orbit 
coupling arises as relativistic quantum effect from the 
Dirac equation, and describes the interaction of the elec- 
tron spin, momentum and electromagnetic field. In the 
system with spin-orbit coupling, the semiclassical equa- 
tions of electron were studied recently, and some novel 
effects have been found [6-15]. However, in the sys- 
tems with spin-orbit coupling, there remain some ques- 
tions to be answered. As well known, the non-relativistic 
approximation of Dirac equation can be obtained from 
Foldy-Wouthuysen transformation. In the transforma- 
tion, there exists a gauge potential in momentum space 
[15,16], so one question is, what is the position operator 
in semiclassical equations? Since the position operator 
is the space-time parameter, its definition is significant. 
If we are uncapable of defining it correctly, we couldn't 
have the complete understanding for spin-orbit coupling. 

In this paper, we investigate the gauge field of posi- 
tion operator in momentum space, and its related effects. 
Based on the dynamic continuity equation, we derive the 
quantum force and torque which contain two parts. The 
conventional part has the same form as in the classic 
electrodynamics, and the topological part originates from 
the spin-orbit coupling and other terms from relativistic 
quantum correction. We notice that the topological part 
has a close relation to the Zitterbewegung effect. For a 
two-dimensional system in a magnetic field, we propose 
that the topological force and torque can reveal more 
complicated phenomena. 

Topological velocity The Dirac equation of electron 
with the wave function = {f,x) T reads ih-^ 1 ^ = Wf>, 
with the Hamiltonian H = col ■ it + /3mc 2 + V(x), where 
7r = P -A, a and j3 are the 4x4 Dirac matrices, 
V(x) = e<j> is a scalar potential, and A is a vector po- 
tential for a magnetic field, B = V x A. m and e 
are the electron mass and charge, and c is the speed of 
light. To reveal the spin-orbit coupling in Dirac equa- 
tion, we perform the Foldy-Wouthuysen transformation 



^ = [/(■7r)\E f , where U(ir) — e , it is a unitary trans- 
formation [17,18]. Choosing S = —i{j3ot ■ Tr/2mc), and 
substituting it into the Dirac equation, we obtain the 
transformed Hamiltonian 

H = f3mc 2 +e lS ca-ne- lS +e lS V(ihd p )e- lS -e lS ihd t e- lS . 

(1) 

The scalar potential becomes V(D) with covariant 
derivative defined by D = ihd p + A, with the pure 
gauge potential A = ihU {izy d v U (li) . By defining the 
covariant derivative, we have the position operator X = 
ihd p + A. The gauge potential A here is trivial. Ne- 
glecting the inter-band transition, and considering the 
adiabatic approximation, we find the non-trivial gauge 
potential 



A = 



h(7T X E) 

4m 2 c 2 



(2) 



where E = 1 ® <x. 

By expanding the Hamiltonian H , and taking S 
hea.-E/Am 2 c 3 to do another transformation 
we get the non-relativist Hamiltonian 
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where the first two terms refer to kinetic energy and 
scalar potential, the third is the Zeeman energy, the 
fourth term is a Darwin term, refers to the correction 
of scalar potential, and the fifth term is so-called the 
spin-orbit coupling term. With this Hamiltonian, the 
Schrodinger equation with wave function ip then can 
be written as ih-j^tp — H sc hip, where tp = (1 + (<x • 
7r) 2 /8m 2 c 2 )</3 is the wave function <p which have been 
performed unitary transformations. After the unitary 
transformations, the wave function x ~ 1/c 3 , so that we 
can neglect it. Based on this fact, we can regard the 
wave function ip have been done the normalization. In 
the Schrodinger equation of ip, the position operator has 
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been redefined by a covariant derivative in the momen- 
tum space X = D = ihd p + A, with the gauge potential 
A = h(-K x cr)/Am 2 c 2 . 

We now turn our attention to the gauge potential in- 
duced by transformation U(tt). In our process, we con- 
sider the adiabatic approximation, so it is straightforward 
to interpret the gauge potential A as the non-Abelian 
Berry gauge potential, and it is Berry connection in mo- 
mentum space, the Berry phase = § c dp- A, where C is 
loop in the momentum space. The Berry gauge potential 
describes the non-trivial geometry of the fiber bundle of 
the Hamiltonian eigenvectors over the phase space, it de- 
scribes the influence of negative energy states in positive 
energy space. 

Using A = h(-7v x er)/4m 2 c 2 , and the relation F^ = 
d Pi Aj — d Pj Ai — (i/h)[Ai, Aj], we get the non-Abelian 
Berry gauge curvature F^ = —Eijkfr(<Jk/2m c — (cr ■ 
7r)7Tfc/877j 4 c 4 ). Let A = — fcfc, tv w P, introducing a 
dual vector Ft — SijkFij, then the gauge field can be 
rewritten as F k = (X/m 2 c 2 )(l - P 2 /4to 2 c 2 ), where the 
first term presents a monopole taking the radial length 
at the rest momentum mc with its strength given by \(Jk 
in the momentum space, the second term is the higher 
order of correction. If we neglect the second term in Fk, 
then 



A 



(mc) 2 



(4) 



We note that this field only has two values ±77 in the 
momentum space, so the difference between this field and 
the ordinary field Fk = A/P 2 [19], is that this field is 
uniform. Obviously, this gauge field is a consequence 
of non-relativistic approximation P <C mc. In classical 
dynamics, the rest energy mc 2 does not affect the motion 
of particle in non-relativistic limit. However, in quantum 
mechanics, the rest momentum mc induces a gauge field 
in the momentum space, then it will influence the motion 
of particle. 

We notice that this gauge field can induce topological 
effects. Noting that after the unitary transformation, the 
covariant position operator of Dirac equation is defined 
by X, and the commutator of covariant operator is not 
trivial: [a;j,a;j] = ihFij. This relation is characteristic 
of the unitary transformation, it is this relation that in- 
duces the gauge field in non-relativistic approximation. 
Defining x = ihd p as the canonical position operator, 
then X = x + A, and the semiclassical equation of mo- 
tion is X = x + tt x T, with T k = Fk/2. Apparently, 
the last term of the equation above is a topological term, 
and results from the non-trivial commutator of covari- 
ant position operator. It acts as the Lorentz force in the 
momentum space. 

By using the Hamiltonian (3), we define the velocity 
X = x+ (eft/4m 2 c 2 )(E+ ^ x B) x cr, where x = -rv/m + 
(eS/4m 2 c 2 )E x cr is the conventional velocity. The second 
term in X actually is a topological term, it is caused by 



the gauge field Fk — A/(mc) 2 , and we regard it as a 
topological velocity. Introducing the electron intrinsic 
magnetic moment m = (eh/2mc)cr, we can rewrite the 
velocity as 



(E x m)/2c 



(5) 



where v = x is the conventional velocity operator which 
corresponds to position operator x, E = E + (jr/mc) x 
B is the total electric field which affects an electron in 
its local coordinate frame. If we remind that the gauge 
field is induced by the interference between the negative 
and positive energy states, this topological term in Eq. 
(5) is contributed by the rest energy under the influence 
of negative states. Comparing with two velocities, we 
find the terms concern with E in Eq. (5) is twice that 
of conventional velocity. Besides, our definition adds a 
contribution of magnetic field. One can easily find that 
these differences between the two velocities is caused by 
the gauge field in momentum space. 

Topological Force and Torque It is very interesting 
in the topological force which acts on the electron in an 
electromagnetic field. Let us start from the dynamic con- 
tinuity equation. The dynamic continuity equation rep- 
resents the change form of a physical quantify in space- 
time, it always requires the existence of a source such as 
a force or a torque. 

Denote by p M — xp^ p M ip the momentum density, 
where p M = mv. Using ^p M tp = ^p M ij; + ^p M tp 
and ihip = H sc hip, then we arrive at 



p M + V • J 



M 



f, 



(6) 



where J M = ^3™^ = [l/2)^{p M , v}4>, and v is 
the conventional velocity. The force is given by f = 
(l/ih)[p M ,H ach ] = i con + Hop, where f co „ is the con- 
ventional force which corresponds to the position oper- 
ator x, and f top is the topological force. The conven- 
tional force can be read f co „ = p^E + (J c x B)/c + 
(VB) • m + (VE) • P + (e/2mc 2 )E x (rh x B'), where 
p c = eip^ip and 3 C — e^vip are the charge and the 
current densities for f co „. m = and P is defined 

by P = ip^Ptjj — ^(2mc x I i 1 )V'! it correlates with the 
electrical polarization induced by the spin-orbit coupling 
[10]. B = B — (7r/2mc) x E presents the total magnetic 
field which an electron experiences in its local coordinate 
frame. Meanwhile, considering the topological velocity, 
the topological force can be found to be 



Hop 



2mc 2 



E 



(rh x B ). 



(7) 



It is clear that the topological force is an additional 
term to the conventional force, so that we can divide the 
total force into two parts: f = f c ; a + f so where f c ; a is the 
force the electron experiencing in classical motion, and 
f so is the force produced by the spin-orbit coupling. The 
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FIG. 1: The electrical field E and the magnetic field B exert 
the pure quantum forces f so on the opposite spin polarized 
electrons, where and jf are charge currents with different 
spin polarization. 

force f c ia can be written as f c ; a = p c E + (J c x B)/c + 
(VB) • m + (VE) • P, it is well known that the classical 
forces are Lorentz, "Stern-Gerlach" , and electric dipolc 
forces. All the forces above have correspondence forces 
in classical dynamics, so we combine them to a classical 
force f c i a . 

The force f so can be read as f so = (e/mc 2 )E x (m x 
B ), where E = E + (iv/2mc) x B. To explore the 
physical meaning, we divide this force into two parts, 

f so = e V^ ^mc™" 1 x B -0 + eip^ - x where the 

first term indicates the electron experiencing an electrical 
field E = [(E" x m)/mc 2 ] x B , and it is an electric 
force. One can find that the second term of f so is a pure 
quantum effect, since (E x B )/mc 2 is dimensionless, this 
term does not concern with any classical quantify, it is a 
pure quantum force, as shown in Fig. 1. 

When the magnetic field vanishes, the force f so reduces 
tof e/e = -i/it(e/2mc 2 )(ih'E)(^xE)f This force origi- 
nates from the second term of f so , so it is a pure quantum 
force, it can be used to explain the Zitterbewegung of the 
electron in spin-orbit coupling system [7,8] . The direction 
of this force is perpendicular to the electrical field, and to 
the velocity and the spin polarization direction, so it is a 
transverse force. When the electrical field vanishes, the 
force f 50 becomes i mag = -^(e/2mc 2 )rh • x B)Bip, 
this force also originates from the second term of f so , 
and it also is a pure quantum force. The direction of this 
force parallels with the magnetic field, and is perpendic- 
ular to the velocity and the spin polarization direction, 
so this force also is a transverse force. In fact, this force 
has duality with f e ; e , so it also can be known from the 
Zitterbewegung effect. 

We now turn to the angular momentum continuity 
equation. Denoting p AM as the total angular momen- 
tum density, then p AM = p° + p s , where the orbit mo- 
mentum density p° — ^X x p M ip, the spin density 
p s = (h/2)ip^crip. Similar to the momentum continuity 
equation, we arrive at 

•p AM + V • J AM = T. (8) 

The angular momentum current is denned by J AM — 
(l/2)-0t{Xx J M + (ft/2)er, v}V>, where the first term is the 
orbit angular momentum current, and the second term 
is the spin current J 5 . The total torque is given by 



T = (l/ih)[p AM , H sch ] = T con + T top , where T con is the 
conventional torque which corresponds to the position 
operator x, T top is the topological torque. The conven- 
tional torque can be read T co „ = ip^x xf-0 + mxB + 
P x E + ip'(TT/2mc) x (E x m)ip, where f is an operator 
with f = tpifip. Meanwhile, considering the topological 
velocity, the topological torque reads 

Ttop - xfip + V^x x f top1 p, (9) 

where the first term concerns with the gauge poten- 
tial, and the second term contributed by the topolog- 
ical force, so both two terms can be known from the 
Zitterbewegung effect. 

Of cause, the total torque can be divided into the or- 
bit torque T° and the spin torque T s . Combining the 
conventional and the topological torque, the orbit torque 
has the form T° = tp^X. x ftp. Obviously, T° has a cor- 
responding torque in classical dynamics. The spin torque 
is given by T s — m x B — tp^m x x E)ip. The first 
term is contributed by the magnetic field. The second 
term is the influence of electric field E, which can in- 
duce a magnetic field in electron's own coordinate frame, 
where 1/2 is Thomas precession factor. These two terms 
can be found in Bargman-Michel-Telcgdi Equation which 
describes the precession of the electron's spin [20]. 

Now we can write the spin continuity equation as p s + 

V • J s — T s , where the spin density is defined by p s = 
(h/tyip^crip, the spin torque is defined by T s . The spin 
current is J s = (h/4)ip^{(T, v}ip, v is the conventional 
velocity operator. 

Nonlinear Spin-Orbit Coupling We apply the pre- 
vious results to a two-dimensional system under a weak 
magnetic field B = (0,0, B), with A = (0,xB,0). 
The Hamiltonian reads H 2 d = n 2 /2m — ub& • B + 
(eF/2)(E X <t) ■ 7r + eEx, where ub is Bohr magnetic 
momentum, F — (h/2)/(mc) 2 is the gauge field, the 
electrical field E = (E, 0, 0). The topological veloc- 
ity v y — hky/m — p,BFBa z k y — eFEa z , where k y is 
wave vector. The spin current operators read J* = 
(h 2 /2m)a z k y - (h/2)fi B FBk y - {h/2)eFE, where the su- 
perscript presents the spin direction, the subscript de- 
notes the direction of electron motion. It is found that 
the topological term here only concerns with momentum. 
In this case, the trajectories read X = x + hFk y cr z /2, 

Y = y. The spin current 

3 z y - a sH E, (10) 

where the spin Hall coefficient <j s h = — (e/4-7r)((5 + (5e/4 — 
e/2), where S = eBF/c, e = hkpF are dimensionless pa- 
rameters, kp is the Fermi momentum. It is clear that the 
second term in the spin Hall coefficient <j s h is contributed 
by the gauge field. 

For spin current J^, there is a pure quantum force 
perpendicular to the magnetic field, it is given by ff = 
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HB{eE /mc 2 )^ (JzipB . We can find that this force is pro- 
portional to the magnetic field. The directions of the 
spin polarization, the current and the force are perpen- 
dicular to each another. This force originates from the 
Zitterbewegung effect. 

If we construct a spin-orbit coupling Hamiltonian as 
follows, 

7T 2 

H = — - /i B <r ■ B + V (tt x 9) ■ u, (11) 
Am 

where r\ = (he/4m 3 c 3 )B is a spin-orbit coupling param- 
eter, v is a unit direction which is along the z direction, 
and 9 = 9(ir,cr) = tt x cr is a vector presenting spin- 
orbit coupling types, we can also obtain the same ve- 
locity of the two-dimensional system above (the system 
has been eliminated the scalar potential) . This indicates 
that the Hamiltonian H2D, which takes the gauge field 
into account, is equivalent to a system with more complex 
spin-orbit coupling. This equivalence illuminates that the 
gauge field can present more complicated spin-orbit cou- 
pling in spin-orbit coupling system, such as tt x 9, with 
9 = tt x er or (7r-er)n, where n is a unit direction. In fact, 
Hamiltonian (11) is a generalization of Rashba model in 
which 9(tz,(t) = cr [21]. 

In the usual spin-orbit coupling system, the types of 
spin-orbit coupling are tt x er and tt ■ cr. They will exist 
as long as the spin-orbit coupling potential is not weak. 
However, these two types only present the first level of 
spin-orbit coupling, the types which we consider present- 
ing the complicated level of spin-orbit coupling, such as 
tt x 9, can also exist in a spin-orbit coupling system. 
Comparing with the first level, the complicated level of 
spin-orbit coupling is the coupling of the orbit and the 
first level. The first level is proportional to the orbit pa- 
rameter 77, and the complicated level is proportional to 
the square of the orbit parameter, so one can find that 
the complicated level reflects the non-linear character of 
spin-orbit coupling system. In fact, there are not only the 
linear spin-orbit coupling, but also the non-linear spin- 
orbit coupling in a spin-orbit coupling system. The first 
level of spin-orbit coupling such as tt x er or tt ■ cr is the 
linear spin-orbit coupling, the complicated level such as 
tt x 9 reflects the non-linear spin-orbit coupling. From 
this point, we say that the complicated level presents the 
non-linear character of a spin-orbit coupling system. 

The significance of the topological force and torque lies 
in the fact that they provide a complete picture of de- 
scribing the motion of a non-relativistic electron in a spin- 
orbit coupling system. It is found that the topological 
force appears as a pure quantum effect when both elec- 
trical and magnetic fields exist in the system, while the 
topological torque originates from the topological force 



and a gauge potential. Through these force and torque, 
the Zitterbewegung effect can be understood deeply, the 
more complicated phenomena such as the nonlinear spin- 
orbit coupling also can be revealed. This implies that the 
spin-orbit coupling can produce more novel effects. 
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